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ABSTRACT 

We present extensive spectroscopy and photometry of the dwarf nova V630 Cas- 
siopeiae. A late-type (K4-5) absorption spectrum is easily detectable, from which we 
derive the orbital parameters. We find a spectroscopic period of P — 2.56387±4 x 10~ 5 
days and a semiamplitude of K% — 132.9 ± 4.0 km s _1 . The resulting mass func- 
tion, which is a firm lower limit on the mass of the white dwarf, is then f(M) — 
0.624 ± 0.056 Mq. The secondary star is a "stripped giant", and using relations be- 
tween the core mass and the luminosity and the core mass and the radius we derive 
a lower limit of M2 > 0.165 Mq for the secondary star. The rotational velocity of the 
secondary star is not resolved in our spectra and we place a limit of Vi-ot sin i < 40 
km s . The long-term light curve shows variations of up to 0.4 mag on short (1-5 
days) time scales, and variations of 0.2-0.4 mag on longer (3-9 months) time scales. In 
spite of these variations, the ellipsoidal light curve of the secondary star is easily seen 
when the data are folded on the spectroscopic ephemeris. Ellipsoidal models fit to the 
mean light curve give an inclination in the range 66.96 < i < 78.08 degrees (90 per 
cent confidence). This inclination range, and the requirement that M2 > 0.165 Mq 
and V T0 t sin i < 40 km s _1 yields a white dwarf mass of Mi = 0.977^0 098 Mq and a 
secondary star mass of Mi = 0.172^Qg 1 2 Mq (90 per cent confidence limits). Our find- 
ings confirm the suggestion of Warner (1994), namely that V630 Cas is rare example 
of a dwarf nova with a long orbital period. 
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1 INTRODUCTION 

V630 Cassiopeiae is a poorly studied cataclysmic variable 
star. This star had a 4 



mag in outburst in 1950 ( Whitney 
197; ) , and a 2 mag outburst observed in late 1992 flHonoy ' 



cutt et al. 1995). The optical spectrum in quiescence shows 



absorption lines indicati ve of a late K-star and l ow-level ex- 
citation emission lines (szkody & Howell 1992). V630 Cas 



is classified in the Duerbeck (1987) catalog as a likely WZ 
Sagittae-type dwarf nova, a particular subtype of outburst- 
ing cataclysmic variable stars which have relatively large 



* Based on observations collected at the William Herschel Tele- 
scope operated on the island of La Palma by the Isaac Newton 
Group in the Spanish Observatorio del Roque de los Muchachos 
of the Institute de Astrofisica de Canarias 

t Visiting Astronomer at Kitt Peak National Observatory 
(KPNO), which is operated by AURA, Inc., under a cooperative 
agreement with the National Science Foundation. 



amplitude outbursts and relatively long intervals between 
outbursts. However, based on its similarity to BV Cen, GK 
Per, and V1017 Sgr, Warner (1994) suggested that V630 Cas 
is a rare example of a dwarf nova with a long orbital period. 
Warner (1994) predicted a period of ~ 6 days, based on a 
correlation between the orbital period and the recurrence 
time between outbursts, and on a correlation between the 
orbital period and the rate of rise to outburst. We have col- 
lected spectroscopic observations of the source over a num- 
ber of years, and from these we show that the orbital period 
is indeed relatively long with P = 2.56387±0.00004 days. We 
also present a long-term V-band light curve of the source. 
Our observations allow us to derive reasonably precise com- 
ponent masses and other interesting binary parameters. We 
discuss our observations and analysis, and the implications 
of these below. 
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Figure 1. The Roboscope V-band light curve of V630 Cas, com- 
plete until November, 2000. 



2 OBSERVATIONS 
2.1 Photometry 

This V630 Cas photomet ry was obtained by RoboScope 
( |Honeycutt fc Turner 1992- ), an unattended 0.41-m telescope 
in central Indiana that is devoted to long-term photometric 
monitoring of accretion systems. The data were red uced us- 



ing the method of incomplete ensemble photometry ( Honey 



cutt 1992[ ). The ensemble solution used 76 comparison stars 
and 579 CCD exposures obtained over the interval 1991 July 
to 2000 November; the zero-point was established using 12 
secondary standards from Henden & Honeycutt (1995). Fig. 
|l| shows the complete light curve, including an extended 3- 
month long outburst in 1992. This 2 mag outburst was de- 
scribed in Honeycutt et al. (1993), and further discussed in 
Warner (1994). We also see in Fig. |l| a rise of ~0.2 mag over 
the four years following the outburst. This is followed by a 
similar -~0.2 mag decline during the succeeding four years. 
Superimposed on this 1993-96 rise and 1997-2000 decline is 
an additional ~0.2 mag of scatter that substantially exceeds 
the observational error. There are changes of up to 0.4 mag 
that are unresolved at our typical data spacing of 1 to 5 
days, superimposed on resolved 0.2-0.4 mag changes with 
time scales of 3 to 9 months. 



2.2 Spectroscopy 

Table ^ is a journal of the spectroscopic observations. The 
KPNO 4m spectra were taken with the R/C spectrograph 
and a CCD detector, and covered AA 4450 - 6930 with ~ 3 A 
resolution; individual exposures were typically 600 s. At the 
MDM Hiltner 2.4m telescope we used the modular spectro- 
graph and a grating covering 4300 - 7400 A with 3 A resolu- 
tion, but with considerable vignetting toward the ends. The 
WIYN spectra were from the Hydra multi- fibre instrument; 
these had inferior signal-to-noise and resolution, but proved 
useful in confirming the ephemeris (below). The spectra from 
the 4.2m William Hershel Telescope (WHT) were obtained 
through the service programme of the Isaac Newton Group 
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Figure 2. Top: The averaged spectrum (in the rest frame of 
the secondary) from the 1994 KPNO observations, normalised to 
unity at 5500 A. Bottom: The spectrum of HD 333388 (a K5 
dwarf) observed with the same instrumentation. The spectrum 
was normalised to unity at 5500 A and scaled by 0.38. Middle: 
The difference spectrum, which is essentially the spectrum of the 
accretion disc. 
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Figure 3. The results of the period search on the absorption line 
velocities. 



of telescopes on La Palma. The instrumental configuration 
consisted of the red arm of ISIS and the R1200R grating, 
yielding « 1 A resolution near Ha. The two spectra each 
consist of the median of three consecutive 20-minute expo- 
sures. We were careful to maintain an accurate wavelength 
calibration for all the observations. We reduced the raw 
CCD pictures to one-dimensional, wavelength-binned, flux- 
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Table 1. Journal of Spectroscopy 
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N is the number of velocities included in the analysis; for the MDM data, each velocity represents the sum of several exposures, 
typically of 8 min. duration. Times given are heliocentric times of mid-exposure of the first and last velocity each night. 



3 ANALYSIS 
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Figure 4. The phased absorption line and emission line velocities 
and the corresponding best-fitting sinusoids. The plot symbols are 
filled squares: KPNO 4-meter absorption line; open circles: WIYN 
Hydra absorption line; filled triangles: MDM 2.4m absorption line; 
5-pointed star: WHT absorption line; small round dots: KPNO 
4-meter H-alpha emission velocities; and small x's: MDM H-alpha 
emission velocities. 



calibrated spectra using standard techniques implemented 
in IRAF0 



3.1 Spectroscopic Period 

Our spectra of V630 Cas show Balmer emission lines, emis- 
sion lines of helium (most notably the lines at A5876 and 
A6678), and relatively strong absorption lines from a «K 
type star (see Fig. ^). Similar features were also noted by 
Szkody & Howell (1992). To measure the absorption spec- 
trum radial velocities, we obtained spectra of K-type stars 
using the same equipment, and cross-correlated thes e with 
the V630 Cas spectra usi ng the IRAF tasks f xcor (Tonry 
Davis L97gb and xcsao ([Kurtz fe Mink 1998|). The cross- 



correlations mostly yielded formal error estimates < 10 km 
s _1 , but the mutual agreement of the standard stars sug- 
gests that systematic effects limited the uncertainty of a 
single velocity to ~ 15 km s _1 . 

We searched for periods in the absorption velocities by 
fitting them with sinusoids over a sufficiently dense grid of 
test frequencies, from 0.01 to over 4 cycle d _1 . Fig. ^ shows 
the result, plotted as the inverse of the mean square residual 
of the fit. This period-finding method should be especially 
suitable for the present data, in which the periodic signal is 
expected to be accurately sinusoidal and the modulation is 
expected to greatly exceed the noise. In this circumstance 
a very large spike is expected at the correct frequency, be- 
cause the small scatter at that frequency results in a much 
tighter fit than at incorrect frequencies. This expectation is 



t IRAF is distributed by the National Optical Astronomy Obser- Research in Astronomy, Inc., under cooperative agreement with 

vatories, which are operated by the Association of University for the National Science Foundation 
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Table 2. Spectroscopic elements 



parameter 



value 



Period, absorption lines (days) 
K2 velocity, absorption lines (km s~ 1 ) 
72 velocity, absorption lines (km s _1 ) 
T , absorption lines (HJD 2,400,000+) 

Period, emission lines (days) 
K\ velocity, emission lines (km s — 1 ) 
71 velocity, emission lines (km s _1 ) 
T , emission lines (HJD 2,400,000+) 



2.56387 + 4 x 10~ 5 
132.9 + 4.0 
-87.3 + 2.3 
49994.4340 ± 0.0087 

2.56375 ± 1.2 x 10~ 4 
39.1 + 4.9 
-111.4 + 2.7 
49662.566 ± 0.034 



borne out here - a frequency near 0.39 cycle d _1 stands out 
prominently, reflecting a uniquely good fit near 2.56 d. The 
single sharp spike shows that the cycle count is uniquely de- 
termined over the 6- year span of the observations. Table [| 
gives parameters of the best sinusoidal fit to the observa- 
tions, and the folded radial velocity curve is shown in Fig. 

The emission lines were generally single-peaked, with 
the exception of the spectra from 1994 November 7, which 
had double-peaked profiles with an average peak-to-peak 
separation of 410+35 km s _1 . We measured the radial veloci- 
ties of the Ha emission lines using the double-Gaussian tech- 
nique of Shafter (1983). The Gaussians had a full width at 
half maximum of 8 A and the separation between them was 
24 A. The emission line radial velocities are shown folded 
on the orbital period in Fig. There is a great deal of scat- 
ter, but the emission line radial velocities are generally out 
of phase with respect to the absorption line radial veloci- 
ties. The formal semi-amplitude of the emission line radial 
velocity curve is K\ = 39 km s _1 . However, we would not 
rely on the emission line radial velocity curve to derive any 
dynamical information given the large scatter seen in the 
velocities. 



3.2 Phase Averaged Light Curve 

We have folded the post-outburst data on the spectroscopic 
ephemeris. Before performing this folding, the post-outburst 
data for each year were normalised to a magnitude of 16.5 to 
remove the long multi-year trend seen in the post-outburst 
data of Fig. |l[ This decreases the scatter in the folded light 
curve but does not correct for the 3-9 month changes seen 
in Fig. |l[ Fig. |s] shows this folded light curve with super- 
posed averages over phase bins 0.05 phase units wide. We 
also constructed a binned light curve by computing the me- 
dian ^-magnitude within each of 20 phase bins using the un- 
normalised light curve (excluding the outburst). The result- 
ing binned light curve is statistically identical to the binned 
light curve shown in Fig. ^| Thus we are reasonably confi- 
dent that the binned light curve shown in Fig. ^| represents 
the true long-term averaged light curve. 



3.3 Spectral type of the Secondary 

We used the technique outlined in Marsh, Robinson, & 
Wood (1994) to find the best fitting spectral type of the 
secondary star and to derive the spectrum of the accretion 
disc. For this we used 26 flux calibrated spectra from 1994 



X3 
=1 

C 

id 







; 










■ ■ I-I 1 V.t" .' , 1 








i , ; 










V630 Cas 






P = 2.56386d 

1 , , 





Phase 

Figure 5. The folded post-outburst V-band light curve (small 
dots). The large filled circles are averages over phase bins 0.05 
phase units wide. 



November, observed with the KPNO 4m. Each spectrum was 
Doppler shifted to zero velocity (using the spectroscopic el- 
ements in Table ^ before the averaged spectrum was made. 
This averaged spectrum, shown as the top spectrum in Fig. 
|j| has a signal-to-noise ratio of » 100 over much of the 
wavelength range. We have 31 observations of 21 different 
stars (mostly dwarfs) observed with the KPNO 4m (13 spec- 
tra were from 1993 and 18 spectra were from 1994). The 
restframe spectrum and the template spectra were all nor- 
malised to unity at 5500 A. Then each template was scaled 
by a weight factor w, and subtracted from the rest frame 
spectrum. The scatter in the difference spectrum is mea- 
sured by computing the rms from a fourth-order polynomial 
fit to the 676 points in the wavelength intervals 5200-5800 A 
and 5950-6500 A. The template spectrum which gives the 
lowest overall rms, corresponding to the smoothest difference 
spectrum, is deemed to be the best match for the spectrum 
of the secondary star. 

The secondary must be evolved well off the main se- 
quence owing to the long orbital period. We have only three 
luminosity class III templates (the rest are luminosity class 
V). However, our experience with this decomposition tech- 
nique indicates that it is not sensitive to the luminosity class 
of the template at these low resolutions. Fig. ^ shows the rms 
values as a function of the spectral type for all of the tem- 
plates. The best two matches were provided by two different 
observations of HD 333388, which classified in the SIMBAD 
database as K8V. However, we believe HD 333388 has been 
misclassified. We adopt K4-5V based on a comparison with 
our other templates, most notably spectra of 61 Cyg A (a 
K5V standard) and 61 Cyg B (a K7V standard). Apart from 
HD 333388, the rms values tend to be the lowest near K4 
and K5. The K7 and later templates clearly give worse fits, 
and the templates earlier than K2 also give relatively poor 
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Figure 6. The rms vs. spectral type for the 31 template spectra 
(see text). The two lowest rms values, and hence the best spectral 
type match, are for the two observations of HD 333388, a K4-5 
dwarf (filled circles). The stars are the three luminosity class III 
templates, and the triangles are the luminosity class V templates. 



fits, although the difference is not as dramatic as for the 
later templates. We adopt K4-5 for the spectral type of the 
secondary in V630 Cas. 

Fig. H shows the results of the decomposition. The top 
spectrum is the observed normalised rest frame spectrum. 
The bottom spectrum is the normalised spectrum of HD 
333388, scaled by 0.38. The spectrum in the middle is the 
difference between the two. The stellar absorption lines are 
removed reasonably well, and the only remaining features 
are the emission lines, the interstellar Na D line, and the 
atmospheric feature near 6870 A. The emission lines near 
4921 and 5015 A seen in the difference spectrum are proba- 
bly due to Hel. The weak emission line near 5169 A might 
be part of the Fe II multiplet 42 (the other two lines would 
be at 4924 and 5018 A, Moore 1972). 

We repeated this general prodecure using the higher 
resolution spectra in an attempt to measure the rotational 
velocity of the secondary star. In this case we have only 
one template spectrum, that of HD 3765, a K2V star. The 
template spectrum was broadened by various trial values of 
V TO t sin i using the standard analytic rotational broadening 
kernel (e.g. Gray 1992) with a linear limb darkening coef- 
ficient of 0.6 and compared to the spectrum of V630 Cas. 
The rotational broadening is not resolved, and we place an 
upper limit of V ro t sini < 40 km s _1 . 



3.4 Light curve models 

The binned light curve shown in Fig. |E] has the character- 
istics of the well-known ellipsoidal variations, namely two 
maxima of roughly equal height per orbital cycle and two 
minima of unequal depth per orbital cycle. The character 
of the ellipsoidal variations is a strong function of the in- 
clination i and the mass ratio Q. Therefore, in principle, 
one can determine the orbital inclination and mass ratio 



by modelling the light curve. In practice, however, there are 
difficulties which can introduce systematic errors into the el- 
lipsoidal modelling. In the case of V630 Cas, the light from 
the accretion disc dilutes the observed ellipsoidal modula- 
tions of the secondary star. The primary minimum in the 
light averaged light curve is relatively deep (it is more than 
0.1 mag deeper than the secondary minimum); this larger 
than expected depth is almost certainly due to a grazing 
eclipse of the star by the disc. There is also a great deal of 
variability that is not associated with the ellipsoidal modu- 
lations. However, with the large number of observations the 
short-term excursions all seem to "average out" , and we are 
reasonably confident that the mean light curve represents 
the true ellipsoidal light curve (which is of course diluted by 
the disc). 



We used the ELC code flOrosz fe Hauschildt 2000[ ) to 
model the binned light curve. We use specific intensities from 
the NextGen models (Hauschildt et al. 1999a,b). The ELC 
code can include light from a flared accretion disc, and geo- 
metrical effects due to the disc (i.e. eclipses of the secondary 
star) are accounted for. We have the usual geometrical pa- 
rameters for the ellipsoidal model, namely the inclination i, 
the mass ratio Q, and the orbital separation a. The main 
parameters for the secondary star are its filling factor /, its 
rotational velocity with respect to synchronous Q, its grav- 
ity darkening exponent f3, its albedo A, and its mean (or 
effective) temperature T tc ff . We assume the star is exactly 
Roche lobe filling and rotating synchronously, hence / = 1 
and Q — 1. We set (3 = 0.08 and A = 0.5, the usual values for 
a star with a convective outer envelope. Finally, the spectral 
type of the secondary star in V630 Cas is K4-5. The tem- 
perature corresponding to this depends on the gravity. For 
example, a K4V star (log g = 4.60) nominally has a temper- 
ature of 4791 K, whereas a K4III star (log g — 1.7) nominally 
has a temperature of 4113 K, according to the parameters 
given in Gray (1992). A simple two-dimensional interpola- 
tion gives T e ff = 4310 K for a K4.5 star with logg = 3 (see 
below). The parameters describing the disc are its inner ra- 
dius -Rinnor, its outer radius -Router, the flaring angle of the 
rim /3 r im, the temperature of the inner edge Idi sc , and the 
power- law exponent of the disc temperature profile £, where 
T(r) = Tdi sc (r/rinncr) ? • Unfortunately, the disc parameters 
are not that well constrained for a number of reasons. For 
example, we have light curves in only one band, so the gen- 
eral slope of the disc spectral energy distribution is rela- 
tively unconstrained. In addition, the bright central parts of 
the disc are not eclipsed, so it is more difficult to pin down 
the disc brightness profile. We therefore fixed the inner and 
outer radii, the flaring angle, and the power law exponent 
at reasonable values (-Rmncr ~ 0.01 Rq, -Router = 75% of the 
primary Roche radius Rrl, /3 r im = 2 degrees, and £ = —0.4, 
a value which gives a flatter profile than the nominal value 
of £ = —3/4 expected for a steady state disc). The only re- 
maining free parameter for the disc in then the temperature 
at its inner edge Tdi sc . Values of Tdi sc near 9000 K give an 
integrated V band intensity roughly equal to the integrated 
V band intensity of the star, which is consistent with the 
observations (e.g. the spectral decomposition shown in Fig. 
|i|). One should bear in mind that the source shows consid- 
erable variability about its mean, hence the "disc fraction" 
measured in Fig. ^| may be slightly larger or smaller than 
the true long-term average. Finally, although the ELC code 
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can include the white dwarf, we omitted the white dwarf 
for simplicity. Since the white dwarf is not eclipsed, its only 
contribution to the light curve would be constant in phase. 
Hence there is no real constraint on its temperature and 
brightness. Since we have only observations in one bandpass, 
we can alter slightly the inner radius and inner temperature 
of the disc to mimic the contribution from the white dwarf. 

We are left with basically four free parameters: the in- 
clination i, the mass ratio Q, the orbital separation a, and 
the temperature of the inner edge of the disc Tdj sc . The two 
binary observables we have are the radial velocity curve of 
the secondary star and the V band light curve. To carry out 
the fitting of the observables, we defined a grid of points in 
the Q-i plane with steps in Q of 0.1 and steps in i of 0.25 
degrees. At each Q,i point in the grid, a and Tdi sc are opti- 
mised using a "grid search" routine adapted from Bevington 
(1969) so that the total \ 2 i s minimised: 



2 

Xtotal 



Xv + Xrv 
20 (v 

( 



mod } i 



60 , 
i=l 



RV m 



The errors on the radial velocities were scaled to yield 
Xrv/ 1 ' = 1 f° r the best-fitting sinusoid. Fig. Q shows the 
contour plot of Xtotal m the Q-i plane. The inclination is 
fairly well constrained to be near 75 degrees, whereas the 
mass ratio Q is basically not constrained (we can put inde- 



pendent constraints on the mass ratio, see Sec. 4.2). There 



is a dramatic and abrupt increase in the \ values near 80 
degrees. At this inclination and higher inclinations substan- 
tial eclipses are predicted, and since none are observed the 
quality of the fits rapidly get worse. At inclinations smaller 
than about 60 degrees the amplitude of the model is far too 
small to match the observations, and the quality of the fits 
also gets much worse. The only observational constraint we 
have from the spectra is that roughly half of the light is not 
from the secondary star. Our models generally satisfy this 
constraint. 

Fig. ^| shows the V light curve and the adopted model 
(solid line) with Q = 5.5 and i = 74 degrees (see Sec. 4.2). 
The fit is good {\ 2 = 21.39 for twenty points). There is a 
grazing eclipse of the star by the outer edge of the disc. The 
dashed line in Fig. |^ shows what the light curve would look 
like if this eclipse is not accounted for (the constant light 
from the disc is still included). The general fitting results 
are not too sensitive on the assumed disc radius. We know 
the disc contributes about half of the light, and if a large 
part of the disc is eclipsed then we would expect a w 0.7 
mag deep eclipse at phase 0. We do not observe such an 
eclipse, which means that the bright central parts of the disc 
are not eclipsed. This in turn constrains the inclination to 
be less than about 80 degrees. In the same vein, the overall 
amplitude of the light curve is relatively large (« 0.25 mag) 
in the presence of substantial disc light. Thus the inclina- 
tion cannot be too small. For most mass ratios, inclinations 
smaller than about 65 degrees produce models with ampli- 
tudes that are too small to match the observations. 




Figure 7. Contours of \ 2 values in the mass ratio-inclination 
plane from the light curve fitting. The filled circle marks the lo- 
cation of the minimum x 2 value. The contours shown correspond 
to the 68.3, 95.4, 99.73, and 99.99 per cent confidence levels for 
two parameters of interest. 




Phase 

Figure 8. The phase-averaged V-band light curve and the ellip- 
soidal model (solid line). The dashed line is a model with light 
from the disc but no eclipses. 



4 ASTROPHYSICAL CONSEQUENCES 

Having laid out all of the observation of V630 Cas, we now 
discuss some of the astrophysical implications of these ob- 
servations. 



4.1 Mass limits on the components 

The observed radial velocity curve of the secondary star K2 
can be used to place a straightforward lower limit on the 
mass of the unseen accreting star M\: 



Mi > /(M) 



PKl Mi sin 3 ; 



2ttG (Mi + M 2 ) 2 
0.624 ± 0.056 M 



(2) 



where M2 is the mass of the secondary star and G is the 
universal constant of gravitation. 
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Figure 9. Contours of the effective temperature T c ff in the core 
mass-envelope plane computed using the relations in Section 



If the secondary star is near the giant branch in the HR 
diagram, we can estimate the mass of its core (and there- 
fore a lower limit on the secondary star M2 ) using three basic 
relations obtained on the evolutionary models of "stripped 
giants". The luminosity and radius of the core are strong 
functions of the core mass M core (e.g. Taam 1983; Web- 
bink, Rappaport, & Savonije 1983; also Phinney & Kulkarni 
1994): 



R 



0.16 M Q 



R 



1.3 







M, 



0.16 M, 



(3) 
(4) 



We also have from Roche geometry and Kepler's Third Law 
a relationship between the mean density of the Roche lobe 
filling secondary star and the orbital period: 



M 2 



3tt 



4/3tt.R| 



P*GRi L (Q)l 



(5) 



where Q = Mi /M2 and Rrl (Q) is the sphere-equivalent ra- 
dius of the Roche lobe for unit separation. For 2 < Q < 20, 
the quantity fl^(Q)(l + ls approximately constant 

(it varies between 9.7 and 10.16). Thus, for our case here, 
the mean density of the secondary star is a function only 
of the orbital period to a good approximation. Therefore 
the Roche radius is a function only of the mass of the sec- 
ondary star. The procedure for finding M corc is simple. For 
an assumed core mass M cola we know its radius R CO ia and 
luminosity L cor c (eq. [],^])- For an assumed envelope mass 
M onv we know the total mass of the secondary M2 and 
hence its radius R2 using eq. (^). The radius of the sec- 
ondary star and the luminosity L corc can then be used to 
find the effective temperature T e «. Fig. ^ shows a contour 
plot of T e g in the M COIC -M cnv plane for situations where 
J?corc < Ri- The contour lines are nearly vertical, meaning 
the core mass M oore is tightly constrained whereas the en- 
velope mass Menv is not. According to Fig. [| T e $ = 4300 
corresponds to M corc > 0.165 Mq. Hence the lower limit on 
the secondary star mass is M 2 > 0.165 Mq. 



4.2 Binary Parameters 

We use a Monte Carlo procedure outlined in Orosz & Wade 
(1999) to derive interesting binary parameters and their un- 
certainties from the results of the light curve fitting. Since 
we have the radial velocity curve of the secondary star, 
each point in the mass ratio-inclination plane corresponds 
uniquely to a set of binary parameters (e.g. component 
masses, orbital separation, etc.). Each point in the mass 
ratio-inclination plane also has a certain probability based 
on the x 2 value of the light curve fit. The Monte Carlo pro- 
cedure selects parameter sets equally from different prob- 
ability intervals and generates probability distributions for 
each interesting binary parameter. Since we don't know the 
mass function (eq. ^|) of the accreting star exactly, 20 values 
of f(M) are drawn randomly from the appropriate Gaussian 
distribution at each Q-i point in the plane. We have three 
additional constraints we add to the Monte Carlo procedure. 
The first constraint, Mi > f(M) = 0.624 M Q , is automati- 
cally satisfied. The second constraint is the one derived from 
the stripped giant evolution, namely M2 > O.165M0. The 
third constraint is the upper limit on the rotational veloc- 
ity, Kotsini < 40 km s _1 . For a Roche lobe filling star in 
synchronous rotation, 



Kotsini _ Rrl(Q) ( 1 



Ki 



(6) 



The upper limit of Vtot sin i < 40 implies Q > 5 for K2 = 
132.9 km s _1 . The results of the Monte Carlo procedure are 
shown in Table ^[ 

The mass of the white dwarf is quite close to 1 Mq , 
while the mass of the secondary star is only slightly above 
the lower limit imposed by the core mass: the hydrogen enve- 
lope has a mass of only « 0.01 Mq. The predicted rotational 
velocity of the secondary star (38 km s _1 ) is just below our 
upper limit of 40 km s _1 . As expected from the relatively 
long orbital period, the orbital separation and secondary 
star radius are relatively large for a CV (8.3 Rq and 2Rq, 
respectively). The small mass and large radius of the sec- 
ondary combine to give a surface gravity of log g — 3.07, 
considerably smaller than the gravities of late-type main se- 
quence K-stars (logg = 4.6, Gray 1992). 



5 DISCUSSION 

The orbital period of 2.56 days is relatively long for a CV. 
The Ritter and Kolb catalog (1998) shows only two CVs 
with longer periods. V1017 Sgr is a dwarf nova with a pe- 
riod of 5.714 days (Sekiguchi 1992), and MR Vel is super- 
soft X-ray source with an orbital period of 3.832903 days 
(Schmidtke et al. 2000). Two other long period CVs are GK 
Per with a period of 1.996803 days (Crampton, Cowley, & 
Fisher 1986) and the old nova X Ser with a period of 1.478 
days (Thorstensen & Taylor 2000). Based on the similarity 
between V630 Cas, GK Per, V1017 Sgr, and also BV Cen, 
Warner (1994) speculated that V630 Cas is a rare example 
of a long-period dwarf nova. Based on correlations between 
the orbital period P and recurrence time T and rate of rise 



logT 
logr 



-0.1 + 1.71ogP 
0.95 + 0.7 log P, 



(7) 
(8) 
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where P is in days, T is in years, and r is in days per magni- 
tude, Warner (1994) predicted an orbital period of w 6 days 
for V630 Cas, more than twice the actual value. A period of 
2.56 days implies r = 17.2 days per mag and T — 3.9 years. 
The observed rate of rise for V630 Cas is about 33 days 
per mag, although this is somewhat poorly defined since the 
two observed outbursts were not very similar. The two out- 
bursts were separated by 42 years, and this is much longer 
than the 3.9 years derived from eq. |^. Indeed, an inspec- 
tion of the RoboScope light curve (Fig. 0) suggests that the 
recurrence time is more than 8 years, although we cannot 
rule out the possibility that a short outburst occurred in a 
gap of the Roboscope coverage (the RoboScope light curve 
shown has 579 observations with the median time between 
observations 1.75 days and the longest time between two ob- 
servations 153 days). It seems that the correlations given by 
eq. |^ have little predictive power at present, mainly because 
they are defined by very few systems at long orbital periods. 



6 SUMMARY 

We have presented extensive photometry and spectroscopy 
of the dwarf nova V630 Cas. The period is relatively long 
at P = 2.56387 ± 4 x 10" 5 days. The observed velocity of 
the secondary star (K 2 = 132.9 ± 4.0 km s _1 ) places a firm 
lower limit of f(M) = 0.624 ± 0.056 Mq for the mass of 
the white dwarf. By modelling the ellipsoidal light curve we 
find an inclination between 67 and 78 degrees (90 per cent 
confidence). Our upper limit on the observed rotational ve- 
locity of the secondary star of V TO tsmi < 40 km s _1 im- 
plies a mass ratio greater than 5. Evolutionary considera- 
tions give a lower limit of 0.165 Mq for the core mass of 
the secondary star. All of these constraints yield component 
masses of Mi = 0.977±^ || Mq for the white dwarf and 
M 2 = 0.1721°' Mq for the secondary star (90 per cent 
confidence limits). 
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Table 3. Astrophysical Parameters for V630 Cas 



parameter 


central value 


90 per cent ranges 


i (degrees) 


74.007 


66.962-78.083 


mass ratio 


5.476 


4 959-6 1 05 


Mi (Mq) 


0.976 


0.880-1.141 


M 2 (Mq) 


0.172 


0.165-0.197 


total mass (Mq) 


1.170 


1.039-1.325 


orbital separation (Rq) 


8.269 


8.024-8.687 


R 2 (Rq) 


2.014 


1.965-2.110 


L 2 (bolometric, Lq) 


1.285 


1.001-1.839 


log 32 (cgs) 


3.071 


3.058-3.090 


Vrot sini (km s ) 


38.464 


37.257-40.000 



Note: Additional constraints are M 2 > 0.165 Mq and V TO t sini < 
40 km s _1 . The quote uncertainties are at 90 per cent confidence. 
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